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Abstract. The puzzle of birth velocities of pulsars (pulsar kicks) could be 
£SJ ■ solved by an asymmetric explosion of type II Supernovae. We propose a sim- 

^ ' pie hydrodynamical mechanism in order to explain this asymmetry, through 

l/") \ the advective-acoustic cycle (Foglizzo 2002) : during the phase of stalled 

shock, an instability based on the cycle between advected perturbations (en- 
tropy / vorticity) and acoustic perturbations can develop between the shock 
and the surface of the nascent neutron star. Eigenfrequencies are computed 
numerically, improving the calculation of Houck & Chevalier (1992). The 
linear instability is dominated by a mode / = 1, as observed in the numerical 
simulations of Blondin et al. (2003) and Scheck et al. (2004). The frequency 
• dependence of the growth rate reveals the presence of the advective-acoustic 

Oh | cycle. 

6 

H 

' 1 Stationary accretion above a solid surface 

We consider a shocked accretion flow onto a solid surface (where the velocity 
• <-h , of the flow is null) at a constant accretion rate. A cooling region is located 
' above the surface and described by the generic function C oc pP~ a P a as in 
^ . Houck & Chevalier (1992), hereafter HC92. The basic equations of the flow are 
the continuity equation, the Euler equation and the entropy equation. This lat- 
ter is : + (v.V)S + j; — 0, where a measure of the entropy is defined by 
S = 1/(7—1) log(P/ p 1 )- These equations are perturbed and projected onto spher- 
ical harmonics YL . 



The jump conditions at the shock r s j, are given by the Rankine-Hugoniot rela- 
tions for the stationary quantities. The perturbations at the shock are evaluated 
taking into account the displacement A£ of the shock position and its velocity 
Av = —iujAt;. In particular, the perturbations of the transverse velocity are as 
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follows (Landau & Lifchitz 1987) : 

V1-V2 5A£ 

Sv 9 = 7^- (1.1) 

- v 2 dA£ 

5v <t> = ^~B~^T L2 ) 

r sh sin a</> 

where t>i and «2 are the pre-shock and post-shock velocities of the flow. 
These transverse velocity perturbations at the shock are not null for non-radial 
perturbations. This contrasts with Eq. (51) of HC92, who did not allow for trans- 
verse velocity perturbations at the shock. 

The eigenfrequency w is a complex number (u> r , u>i) such that the velocity 
perturbation satisfies a wall type condition (Sv/v = 0) at the surface of the ac- 
cretor. The imaginary part uii of the eigenfrequency is the growth rate of the 
perturbations. 

2 Calculations of the eigenmodes 

We performed several calculations of the fundamental modes (an example with 
7 = 5/3, a = 1/2 and (3 — 2 is shown in Fig. 2J. The mode I = 1 is always the 
most unstable. This result differs from the analysis made by HC92 because of the 
error in their boundary conditions at the shock. 

The advective-acoustic instability is based on the cycle between advected per- 
turbations (entropy / vorticity) and acoustic perturbations between the shock and 
the surface. A reference timescale r for this mechanism is equal to the accretion 
time from the shock to the coupling region near the surface plus the time for an 
acoustic wave to reach the shock : 

1 1 dr 

— (2.1) 

l-M\v\ [ ' 

The acoustic time £ ac is defined by the time needed for an acoustic wave to prop- 
agate from the shock to the accretor and then back up to the shock, cj ac = 2-7r/f ac 
being the pulsation associated to this acoustic time : 

rr S h 9 dr 

^ Sr. 

On Figs. HIS the growth rate u>i is at best comparable to o; s h ~ 1/t and the 
pulsation uo r of the fundamental unstable modes is close to 27r/r, as expected in 
the advective-acoustic mechanism. The frequency dependence of the growth rate 
of the eigenmodes shows an oscillatory behaviour with a period comparable to 
w ac (Fig. |5J . Such oscillations are expected in the advective-acoustic instability, 
as a consequence of the modulation of the advective-acoustic cycle by the purely 
acoustic cycle (Foglizzo 2002). This interpretation is confirmed by measuring the 
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ratio T/t ac ~ 4.5 — 6.5 which is comparable to the number of modes found per 
oscillations (Foglizzo 2002). 

We note that for big cavities, the most unstable eigenmodes correspond to low 
frequencies, in the "pseudo-sound" regime (uu r < uj ac ). Unstable eigenmodes are 
also found in the acoustic regime (ui r > w ac ), with a smaller growth rate however 
(Fig.©. 




Xso/r^ Xso/r^ 



Fig. 1. For 7 = 5/3, a = 1/2 and j3 = 2, numerical calculations of the frequency u) T and 
the growth rate uoi in units of u s ^ = — u s h/('"sh — »"*) of the fundamental modes I = 0, 1, 
2, 3 depending on the size of the cavity x s o/r* = (r s h — r*)/r*, as in Houck & Chevalier 
(1992). The degree I of the modes is indicated on each curve. 




Fig. 2. For 7 = 4/3, f3 = 2, a = 1/2, and a cavity r s h/r* ~ 12, numerical calculations 
of the frequency u r in units of oj ac and the growth rate Wj in units of w s h of radial / = 
and non-radial I = 1, 2 eigenmodes. 
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On Fig. |3J for 7 = 4/3 and a cooling function described by a = 6, /3 = 1, rele- 
vant to the phase of a stalled shock in type II Supernovae (Bethe & Wilson 1985), 
an unstable mode / = 1 can also be found for big enough cavities (r s h/r* > 3.5). A 
frequency dependence study, still in progress, shows unstable modes I = 1 both in 
the "pseudo-sound" regime and in the acoustic regime. The radial modes (I = 0) 
are always stable. 
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Fig. 3. For 7 = 4/3, a = 6 and (3 = 1, numerical calculations of the frequency u> r in 
units of oj ac and the growth rate LUi in units of u> s h = —v s h/(r s h — r*), of the first unstable 
mode 1 = 1, depending on the size of the cavity x s o/r* = (r s h — r*)/r«. The vertical 
dashed line correspond to r s h/r* = 3.5 (x s o/r, = 2.5). 



3 Conclusion 

The stalled accretion shock above a neutron star is unstable, with a domination of 
a mode I — 1 if the shock radius is large enough. The instability is interpreted as 
an advective-acoustic cycle modulated by a purely acoustic cycle. The instability 
is also found when the cooling function mimics neutrino cooling in type II Super- 
novae. The advective-acoustic cycle is thus a good candidate to seed an asymmetric 
explosion which could lead to an important birth velocity of the neutron star. 
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